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Hidden charmonium decays of Z c in intermediate meson loops model 

Gang LJ*| 

Department of Physics, Qufu Normal University, Qufu 273165, China 
(Dated: April 17, 2013) 

The BESIII collaboration reported an observation of a charged charmonium-like structure Z c (3900) in the 
invariant mass spectrum of //\|/7t in e + e~ — > J/\\m + Tl~ at the center of mass 4.260 MeV, which could be an 
analogue of Z/,(10610) claimed by the Belle Collaboration. In this work, we investigate the hidden charmonium 
transitions of Z c via intermediate D D* meson loops. Our results show that the intermediate D D* meson loops 
are crucial for driving the transitions of Z c . We also make a predication for the hidden charmonium decays of 
Z c 's partner (Z' c ). 

PACS numbers: 13.25.Gv,14.40.Pq,13.75.Lb 



I. INTRODUCTION 



Recently, a new charged state Zf is observed in the J /yi^ invariant mass spectrum of 7(4260) — > J/\\tn + K 
decay by the BESIII Collaboration (H. The reported mass and width are M z ± = 3899.0 ± 3.6 ±4. 9MeV and 
r z ± — 46 ± 1 1 .3 ± 12.6MeV [l|]. Belle Collaboration also observed a new charged charmonium-like structure in 
the ^/\|rTC =t= invariant mass spectrum with 5.2a significance, with mass M z ± = 3894.5 ±6.6±4.5MeV and width 
r z ± = 63 ± 24 ± 26MeV yy. Differing from the other hidden charmonium-like states, such as X(3872), 7(4260) 
etc, Z c is electric charged states. Such a state, if it exists, need at least four quarks as minimal constituents, which 
qq , makes them ideal candidates for exotic hadrons beyond the conventional qq mesons. On the one hand, this chain 

decay mode reminds us of the observations of Z 6 (10610) and Z^(10650) in T(5S) ->■ Z b (Z' b )% -> Y(nS)%K HH. 
On the other hand, the mass of Z c (3900) is in the vicinity of DD* + c.c. This similar phenomenon (mechanism) 
shows that Z L (3900) may be an analogue of Z/,(10610) claimed by the Belle Collaboration. We also expect the 
experimental observation of Z' c and the neutral partners. 

Recently, many investigations have been carried out to explain this exotic states J5l 4l2ll . The results of Ref. 15|] 



show that it is necessary to include explicit Z c (3900) poles, i.e. a resonance structure as an isovector partner of 
X(3872) in order for a more detailed description of the data. Applying heavy quark spin symmetry and heavy 
flavour symmetry the authors found a promising isovector 1 + ~ DD* virtual state near threshold that might very 
k>( ■ well be identified with the newly discovered Z c (3900) lll|2|] 

The intermediate meson loop (IML) transitions, or known as final state interactions (FSI), have been one of the 
important non-perturbative transition mechanisms in many processes II 13142711 . In the energy region of charmonium 
masses, with more and more data from Belle, BaBar, CLEO and BES-III, it is widely recognized that intermediate 
hadron loops may be closely related to a lot of non-perturbative phenomena observed in experiment 11161 - 13 Oil . 



e.g. apparent OZI-rule violations, sizeable non-DD decay branching ratios for \|/(3770), and the HSR violations in 
charmonium decays. Recently, the IML transitions are also applied to bottomnium decays J3ll435ll . By applying the 
on-shell approximation, the bottom meson loops were suggested to play an important role in the Y(5S) transitions 
to the lower T states with the emission of two pions 13111 or one r\ 113211 . This mechanism seems to explain many 
unusual properties that make the T(5S) different from T(45). Similar approach was also applied to the study of 
Zj, and Z' b by Liu et al. 113311 . Within a nonrelativistic effective field theory, the decays of the Z/,(10610) and the 



Z/,(10650) to T(tiS)k and lit,(mP)% are investigated in Ref. [34-] . The power counting analysis in Ref. [34] shows 
that the triangle transition Zj, — > hi,(mP)n is not suppressed compared to Zj, — > Y(nS)n, although the decay is via 
a P-wave. In Ref. J35I1 . we investigated the transitions from the Zi,(10610) and Zi(10650) to bottomonium states 
with emission of a pion via intermediate BB* meson loops in the ELA. By investigating the a-dependence of 
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partial decay widths and ratios between different decay channels, we show that the intermediate BB* meson loops 
are crucial for driving the transitions of Zi,/Z' b — > T(nS)n with n = 1,2,3, and hb(mP)n with m = 1 and 2. 

Since the Z t (3900) is very close to the DD* threshold, the IML should be an possible mechanism in its decays. 
In this work, we will investigate the decays of Z c ' —> J/\\m, \[/'7t and h c % via intermediate D-meson loops in an 
effective Lagrangian approach (ELA) with quantum numbers I G (J PC ) = l + (l + ~) for the Z c /Z' c . The paper is 
organized as follows. In SecHU we will introduce the formulae for the ELA. In Sec. [TTTl the numerical results are 
presented. The Summary will be given in Sec. HVl 





FIG. 1: The hadron-level diagrams for Z+ — ► \|nt + with d(*>D(*> as the intermediate states. Similar diagrams for Z c and Z® 
states decays. 




FIG. 2: The hadron-level diagrams for Z+ — ► h c Jt + with D^*'D^*' as the intermediate states. Similar diagrams for Z c and Zj? 
states decays. 



II. TRANSITION AMPLITUDE 



In order to calculate the leading contributions from the bottomed meson loops, we need the lea ding order effec- 
tive Lagrangians for the couplings. Based on the heavy quark symmetry and chiral symmetry IJ36L 13711 . the relevant 
effective Lagrangians used in this work are as follows, 

V W ^ W = igyDDV^DD-Dd^-g^DZ^d^vidaDlD + Ddabl) 



4w>MdM = gh^DK{DDl + D;p) + ig hcD , D ^^d,h cv D* a Dl , 

t^DM* = -igirmiDid^PijDf -D*%PijDj) + ^ir^^D^Pij V a D* p , 

Li,) M M = g w zi' )fi (DD* u + D* u D)+ig („ e^ va %Z { c (JD* a D*o , 

2£>d(*)d(*> °zVdd*d k f f ' °z { J'd*d* f cv a P' 



(2) 
(3) 

(4) 



where DM = (d^ + ,D^°J and D^ T = LdM-,J5M°J correspond to the charmed meson isodoublets. The 
coupling constants will be determined in the following. 

The loop transition amplitudes for the transitions in Figs. Q] and [2] can be expressed in a general form in the 
effective Lagrangian approach as follows: 



M f 



d 4 q 2 v TiT 2 T iT1 2 

271 D* pol fll02fl3 ' 



(5) 



where T/ (i = 1,2,3) are the vertex functions; a,- = qj — mf (i = 1,2,3) are the denominators of the intermediate 
meson propagators. We adopt the form factor, Ylf!fi(mi,qf), which is a product of monopole form factors for each 
internal mesons, i.e. 

]j!Fi(mi,qj) = Tx(mi,q\)<F % (rm,q\)fr(m^q\) , (6) 

i 

with 

2 , /A?-/n?\ 
f i{mi ,qt)={^—ty (7) 

where A, = m, + ocAqcd and the QCD energy scale Aqcd = 220 MeV. This form factor is supposed to parameterize 
the non-local effects of the vertex functions and remove the loop integral divergence. 

The explicit transition amplitudes for Z c (pi) -> D^(q\)D^(q3)[D^*\q2)] -> V(pf)K(px) via those triangle 
loops are given as follows: 

M DD *[ D ] = (if J j^^[gZ b D*D£i fl }[gyDD£* f P {qi-q2)p}[gD*DKPiie} 

— 2~ 2 2 — ~T \\m m i,qi) (8) 

q\ — m\ q^ — m\ q^ — m^ ,■ 

[-so- 

jFi(mi,qf) (9) 



M DD * [D *] = (if J J^^[gZ c D*D^[g^D*D^pc^P P f^fq2\[-gD*D*T l ^KXPlq\] 

i(-g x » + qlqVmi)i(-g^ + <fA/mi) 
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Md*D[D*] = (if I —-^ [gZ c D*DEin] [gM/D*D* (gpag^x - gpigoi; + gp^gcz)^*/ (q\ + <?2) X ] [-gB*B%P%$\ 



(2n) 



q j — m j q^ — m^ q\ — m| ,■ 



M D « D *\ D] = (i) 
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r J 4 

M D * D .p*} = (if J j^^[gZ c D*D*£ M vap>P^^}[gyD*D*(gpag^~ gpxgo^+ gp^gaz)z* f P (qi + q2f][-gD*D*%^KKPlq\] 

f(-g^ + g ? g fM) .■(-g° e + ^/ m j) ii-gto + fiqj/m*) ^ 2 

x 2 2 2 2 2 2 ~~ llmmi,qi) ■ (12) 

(/■J -WJ| (^ — m 2 <?3 — m 3 i 

Also, the explicit transition amplitudes for Z c (pi) ->Z)W(?i)Z)W(j3)[DW({2)] — > h c (pf)n(p n ) via those triangle 





FIG. 3 : (a) The a-dependence of the partial widths of Z+ — » J/\\rn + (solid line) and \|/'7t + (dashed line), (b) The cc-dependence 
of the partial widths of Z' c + — » J/\\m + (solid line), \|/7t + (dashed line). The mass of Z' c is fixed to be 4040 MeV. 



loops are given as follows: 
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(13) 



(14) 



(15) 



Y\7i(m„qf) 7 (16) 



where /?,, p/, p% are the four-vector momenta of the initial Z c , final state charmonium and pion, respectively, and 
q\, q2, and qj, are the four- vector momenta of the intermediate charmed mesons as defined in Figs.Q]and|2] 

III. RESULTS 

Before proceeding to the numerical results, we first discuss the parameters, such as the coupling constants, in 
the formulation given in Section. [TT] In Eq. (fTJ, the following couplings are adopted in the numerical calculations, 



g\/DD = 2g2y/m~ym D , g\s/D*D 



_gj£D_ 

y/m D m D - 



g\s/D*D* = gM/D*D 



m D * 
m D 



-m D * 



(17) 



In principle, couplings g 2 must be computed by nonperturbative methods. It shows that vector meson dominance 
(VMD) would provide an estimate of these quantities [36J]. The coupling g 2 can be related to the J/\\i leptonic 





FIG. 4: (a) The a-dependence of the partial widths of Z+ — > h c n + . (b) The (X-dependence of the partial widths of Z' c + — s> h c n^ 
The mass of Z' is fixed to be 4040 MeV. 
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FIG. 5: (a) The partial widths of Z+ — > y/\|/7l + (solid line) and \|/7t + (dashed line) in terms of the mass of Z' c . (b) The partial 
widths of Z', + — > h c n + in terms of the mass of Z'.. The cutoff parameter a = 0.8. 



constant /y which is defined by the matrix element (0|cy /i c|.//\|/(/7,e)) = /ymy&", and g2 = y/fn^/^mofyf, where 
/y = 405 ± 14 MeV, and we have applied the relation gyoo = '«y//y- 
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FIG. 6: The a-dependence of the ratios Rz (solid line) and R^' (dashed line) defined in Eq. ( I22t , 



The ratio of the coupling constants gy'DD to gv/DD is fixed as Ref. IJ25I1 

gy'DD 



g\s/DD 



0.9. 



In addition, the coupling constants in Eq. (0 are determined as 

gh c DD* = -2giy/mh c m D m D *, gh c D*D* =2g\ 



m D « 



(18) 



(19) 



with g\ = —^/my^/3/fy^, where m Xrf) and f Xc0 are the mass and decay constant of % c o(lP), respectively 13811 . i.e. 
/ Xr0 = 175±55MeVjM. 

In analogy to what is known about the Zf, states we assume that the width of the Z £ (3900) is saturated by the 
DD* intermediate state. We obtain the coupling constant gz c D*D = 6.85 GeV. The relations 



gZ'D*D=gZ c D*D, g 7 



are applied to extract the couplings for g ih 



'N 



m D m- {l) 

A' 



(20) 



The coupling constants relevant to the pion interactions in Eq. (@]i are 



2,? 



gD'DK = —\Jmumj)* , gD*D*n = 



gD*D% 



(21) 



where f % = 132 MeV is the pion decay constant, the parameters gy respects the relation gv = m p /f % 13711 . We take 



7-1 



X = 0.56GeV-' and g = 0.59 A41H . 

In Fig.[3ja), we plot the a-dependence of the decay partial widths of Z c — > J/\\tK (solid line) and \\t'n (dashed 
line), respectively. A predominant feature is that the a dependence of the partial widths are quite stable, which 
indicates a reasonable cutoff of the ultraviolet contributions by the empirical form factors. As shown in this figure, 
at the same a, the intermediate D-meson loop effects turn out to be more important in Z c — > \\i'k than in Z c — > J/\\ik. 
As a result, a smaller value of a is favored in Z c — > \|/'ti. This is understandable since the mass of \|/' is closer to 
the thresholds of DD* or D*D* than 7/\|/ y42||. Thus, it gives rise to important threshold effects in Z c — >• \y'n + . 

One also notices that the a-dependence of the partial widths for Z+ — > \|/7I + are stabler than that for J/\\ik. This 
indicates that the enhanced partial widths are not from the off-shell part of the loop integrals. As we know that the 
dispersive contributions become rather model-dependent near threshold, the enhanced (but rather stable in terms 
of a) partial widths for Z+ — » \\t'n + suggest that more stringent dynamic constraints are presumably needed to 
describe the near-threshold phenomena where the local quark-hadron duality has been apparently violated. What 
makes this process different from e.g. \\f' — > h c %° in Ref. 12611 is that there is no cancelations between the charged 
and neutral meson loops. As a consequence, the subleading terms in Refs. 1261 12711 become actually leading 
contributions. 

For the Z' c decays, we carried out the following two analysis to understand this exotic states due to the lack of 
experimental information for Z' c . First, we fix the Z' c mass and then study the cut-off parameter a dependence of 
Z' —> J/\\ik, \\i'k. Second, we fix the a values (Here we choose a = 0.8) and then study the Z'. mass dependence 
of these decays. 

The measured central mass of Z e is 23.7MeV above the D* + D° threshold and 22.2MeV above the D* + D° 
threshold 14211 . So we put the Z[. resonance by the same amount above the D* + D*° threshold which gives its 
expected mass at approximately 4040MeV 14211 . In Fig. 0b), we plot the a-dependence of the decay partial widths 
of Z' c — > J/\\tK (solid line) and \\t'n (dashed line), respectively. For the Z' c —t J/\\tK and \\t'n, the situation becomes 
more complicated. The thresholds of intermediate meson, i.e., m D ( t ) +m fl («) +2aAgco, are very close to the initial 
Z' c mass at small a region, which leads to much more dramatic dependence behavior. The behaviors become more 
stabler at relatively larger a region. It indicates that for sufficiently heavy decaying charmonia and small values of 
a these singularities are located nearby or even inside the physical regime and additionally introduce unphysical 
contributions. So a relatively larger a should be considered in order to avoid introduce unphysical contributions. 



The a-dependence of the partial widths of Z c — > h c K are shown in Fig.|4ja). From this figure, we can see that the 
intermediate meson loop contributions are important in Z c — > h c Tt, Again, the behavior of Z c — > h c % is quite stable, 
which indicates the dominant mechanism driven by the intermediate meson loops. We make an prediction for the 
partial width of Z'. — > h c n, the results are plotted in Fig. Qb) with m z i = 4040 MeV. Again, Z'. — > h c K shows a 
similar behavior to that of Z' c — > J/\\ik and \\i'k. 

For the Z'. decay, we fix the a = 0.8 and vary the mass of Z' c . The predicated results are present in Fig. As 
shown in this diagram, due to the thresholds effects, there is a peak near 4050 MeV. 

It would be interesting to further clarify the uncertainties arising from the introduction of form factors by study- 
ing the a dependence of the ratios between different partial decay widths. For the decays of Z+ /Z' c + — > J/\\tn + , 
we define the following ratios to the partial decay widths of Z+ /Z' c + — > \\t'n + : 

_ T(Z+^J/ W +) T(Z'+^J/ W +) 

Rz ' = r(z+^ ¥ ';t+) ' Rz '- r(z'+^ ¥ ':r+) ' (22) 

which are plotted in Fig. [6] The stabilities of the ratios in terms of a indicate a reasonably controlled cutoff for each 
channels by the form factor. Since the first coupling vertices are the same for those decay channels when taking 
the ratio, the stability of the ratios suggests that the transitions of Z c /Z' c —> J/\\tn and \\i'k are largely driven by the 
open threshold effects via the intermediate D meson loops. Here, we should mention that since the mass of \|/' is 
closer to D^*'D^*' thresholds than J/\\i, so the ratio only reflects the open threshold effects via the intermediate D 
meson loops. The future experimental measurements of Z c /Z' c — >• 7/\|/7i and \|/7t can help us investigate this issue 
deeply. 

In order to understand this, the following analysis is carried out. First, one notices that we have adopted the 
couplings for the h c and \\t to DD* or D*D* in the heavy quark approximation. Since the physical masses for D 
and D* are adopted in the loop integrals, the form factor will introduce unphysical pole contributions of which the 
interferences with the nearby physical poles would lead to model-dependent uncertainties. By assuming Mo* = 
Mo = 1869 MeV and Mo* = Mo = 2010 MeV, namely, by partially restoring the local quark-hadron duality, we 
calculate the partial widths of Zf and Z' c + —> J/\\ik, \\t'n and h c n. We expect that the partial restoration of the local 
quark-hadron duality will significantly lower the partial widths since there will be only one physical pole in the 
loop and the unphysical one can be easily isolated away from the physical one. This is a rather direct demonstration 
of the sensitivity of the meson loop behaviors when close to open threshold and when the dispersive part becomes 
dominant. It shows a relatively stable behavior, especially Z' c decays at the small a region. 

IV. SUMMARY 

In this work, we investigate hidden-charm decays of the newly discovered resonances Z+ via intermediate D- 
meson loops. In this calculation, the quantum numbers of the neutral partners of these two resonances are fixed 
to be I G (J PC ) = l + (l + ~), which has the same favored quantum number of Zl. For Z c decays, our results show 
that the a dependence of the partial widths are quite stable, which indicate the dominant mechanism driven by the 
intermediate meson loops with a fairly well control of the UV contributions. For the predicated Z'. decays, there 
shows dramatic behavior at small a region, however, the behavior becomes quite stable at relatively lager a region. 
Our results show that it will become sensitive to the meson loop contributions when the final state mass threshold 
are close to the intermediate meson thresholds. Namely, the effects from the unphysical pole introduced by the form 
factors would interfere with the nearby physical poles from the internal propagators and lead to model-dependent 
uncertainties. It is also a consequence of the violation of local quark-hadron duality. Such a phenomenon has been 
discussed in Ref. B27I1 . We expect the experiments to search for the Z' c resonance and the decays of Z c — > DD* and 
Z' c — > D*D*, which will help us investigate the Z c decays deeply. 
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